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Abstract 

The scaling behavior of the closed trajectories of a moving particle gener- 
ated by randomly placed rotators or mirrors on a square or triangular lattice 
in the critical region are investigated. We study numerically two scaling func- 
tions: f{x) related to the trajectory length distribution ns and h{x) related 
to the trajectory size Rs (gyration radius) as introduced by Stauffer for the 
percolation problem, where S is the length of a closed trajectory. The scaling 
function f{x) is in most cases found to be symmetric double Gaussians with 
the same characteristic size exponent a = 0.43 ~ 3/7 as was found at critical- 
ity. In contrast to previous assumptions of an exponential dependence of ns 
on S, the Gaussian functions lead to a stretched exponential dependence of ns 
on 5*, ns ~ e"'^'^'^^. However, for the rotator model on the partially occupied 
square lattice, an alternative scaling function near criticality is found, leading 
to a new exponent a' = 1.6 ± 0.3 and a super exponential dependence of ns on 
S. The appearance of the same exponent a = 3/7 describing the behavior at 
and near the critical point is discussed. Our numerical simulations show that 
h{x) is essentially a constant, which depends on the type of lattice and on the 
concentration of the scatterers. 



Key words: Lattice, particle trajectories, percolation , scaling function, critical 
exponents. 
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1 Introduction 



In a previous paper^^^ we studied numerically the scaling behavior of extented 
(possibly infinite) particle trajectories generated by randomly distributed (ro- 
tator or mirror) scatterers on a square or a triangular lattice. We restricted 
ourselves there to the scaling behavior of structural properties of particle tra- 
jectories — such as of the asymptotic number of right rotators on a particle 
trajectory or of the asymptotic number of sites visited by the moving particle 
on its trajectory with different frequencies — that occurred strictly at critical- 
ity. In this paper we will discuss the scaling behavior found in the region near 
criticality, i.e. we study the critical behavior when one approaches the critical 
point from the outside. This involves, in addition to the scaling exponents r, 
df and a, which already appeared in previous papers*^^"^^ , the determination 
of a scaling function defined by: 



where ns is the probability to find a closed trajectory of length S, Cr and 
Cr^ are the concentration (i.e. the fraction) of right scatterers and the critical 
concentration of right scatterers on the lattice, respectively. For both the 
square and the triangular lattices fully occupied by rotators, one can show , 
by mapping the rotator model onto a percolation problem^^'^^ that has been 
solved exactly before(^), that Cr^ = 1/2, t = 15/7, df = 7/4 and a = 3/7. The 
scaling function f{x) yields more detailed information about the trajectory 
size distribution than the exponent a, where x = {Cr — CrJS'^ . Since there 
are no infinitely extended trajectories away from criticality, f{x) must vanish 
when X — > oo. The exponent — r -|- 1 instead of the usual r occurs because 
the trajectories are constructed here from, what corresponds in percolation 
theory to, a seed and the number of choices to put the seed on a trajectory 
is proportional to the size of the trajectory itself'-^'^'^'^''-*. When x = and 



^5 = S-^+'f[{CR - CrJS''] 
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f{x) becomes a constant, we have, from Eq. (|l]) , 77,5 ~ 5*^^+^, recovering the 
results obtained previouly at criticahty'-^'^-' . 

One important feature about the scahng function f{x) is that it is essen- 
tially invariant once we get into the critical region, i.e. the scaling function 
we obtain at one concentration near criticality is essentially the same as that 
obtained at another concentration. Note that since x = {Cr — Cr^)S" , to 
make f{x) invariant, i.e. concentration independent, we need to give Cr^ and 
a each a unique value. This provides an alternative method to compute both 
the exponent a and the critical concentration Cr^, just from numerical data 
obtained at a few different concentrations in the critical region. 

In addition to /(a;), there is another scaling function h{x) introduced by 
Stauffer'^-'^^^ some time ago. 



Here Rs is the gyration radius of a closed trajectory of length S, in other 
words, -R| is defined as the mean square displacement from the origin of all 
trajectories of length S. The scaling function h{x) characterizes corrections to 
the fractal dimension as one moves away from criticality. From our numerical 
simulations, we found, however, that h{x) is essentially a constant for all x. 
This implies that Rs ^ S'^^^f holds quite accurately in practice even for small 
trajectories, i.e. away from criticality. 

Previous numerical determinations of a were based on the calculation of 
the first moment of the trajectory length distribution <5'>, which is 
divergent at criticality. 
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Since < S >~ {Cr — Cr^)^"' by definition, one has 7 = (3 — r)/cr = 2. In 
this paper we compute the average of the mean square displacement of all 

trajectories, <R^>, which also diverges at criticality, 

00 

<R^> = Y^^l^s 
s 

00 

= Y.S-'^'^^"> f{{CR - CnJSnh'iiCR - CnJSn 
s 

~ {Cn-CRj-^'^'/'f-^^'^ (4) 

If we define an exponent phj <R^ >~ {Cr — CrJ^^ one has p = {2 + 2/df — 
T)jo = 7/3 = 2.333. The hyperscaling relation in Eq. is independent of the 
specific form of j[x) and h{x). Note that the upper limit of the summations 
in Eq. and Eq. (^) is infinity, which poses a difficulty in getting very close 
to the critical point, since the dominant closed trajectories are then too large 
to generate on the computer. On the other hand, we can still compute the 
scaling function j[x) by making a large cutoff (2^^) in S", so that trajectories 
whose length is bigger than 2^^ are not considered. Although then the tail of 
f{x) will get truncated by the cutoff, the remaining part still yields sufficient 
information to determine a . 

The numerical algorithm used to generate the particle trajectories here 
has been explained in the previous paper'^^\ This a Igorithm is very powerful 
due to its speed and its ability to generate large particle trajectories by using 
a virtual lattice scheme and a dynamic memory allocation technique. The 
calculation of each scaling function ][x) or h[x) and the calculation of < S'> or 
<R^> at each concentration involve 500, 000 and 30, 000 particle trajectories, 
respectively. 

The plan of this paper is the following. In section 2 the critical region 
for trajectories generated by rotators on a square lattice is discussed. For the 
critical behavior near Cr^ = Cl^ = 1/2, two qualitatively different scaling 
functions, f{x) a double Gaussian and f'{x) an exponential function, corre- 
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spending to two different values of a: a = 0.43 ~ 3/7 and a' = 1.6, re- 
spectively, near this critical point have been determined, depending on how 
one approaches this point in the phase diagram. While a = 0.43 leads to a 
stretched exponential behavior of ns along the line Ch + Cl = 1, cr' = 1.6 leads 
to, what one could call, a super-exponential behavior along the line Cr = Cl- 
h{x) is obtained by computing f{x)h{x), which appears to be proportional to 
f{x), so that h{x) is essentially a constant. In section 3 we discuss the mirror 
model on the square lattice and a so-called quasi-rotator model deduced from 
the mirror model. In section 4, the critical region on the triangular lattice is 
discussed, which obtains both for rotators and mirrors. We found that the 
exponent a along the critical line is the same as that for C = 1, i.e. a = 0.43, 
and that the scaling functions f{x) and h{x) are qualitatively the same and 
differ only by the values of the constants occurring in them. A discussion of 
the results obtained in this and the previous paper is given in section 5. 

2 Critical region for the rotator model on the 
square lattice 

2.1 The fully occupied lattice 

As we have shown in the previous paper for C = 1, the trajectories generated 
by a moving particle through randomly placed rotators can be mapped onto 
the perimeters of bond percolation clusters. Therefore the theory that has 
been developed for this percolation problem can be applied directly here. The 
critical concentration is C/?^ = C^^ = 1/2 and the exact values for a, dj 
and T are 3/7, 7/4 and 15/7, respectively. From these known exponents and 
the critical concentration, we can compute the scaling function directly from 



Eq. ®, 




S-r+l 
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= S^/'ns (5) 

Note that since ns is the probabihty to find a trajectory of length S, the right 
hand side of Eq. (j^) is an average of S^^"^ , taken over trajectories of length S 
and can be easily determined in our numerical simulations. It is not difficult 
to see that f{x) must be symmetric with respect to x = 0, since the proba- 
bility to generate the same trajectory is invariant under the transformation of 
interchanging Cr and Cl, i.e. 



mcn - \)S'^) = mc, - ^)sn = - ^)sn (6) 

Our numerical calculations of the scaling function f{x) were carried out at 
Cr = 0.47, 0.48 and 0.49, respectively. We found that the scahng functions 
obtained at these three concentrations collapse to a single curve, which could 
be fitted to a double Gaussian, i.e a sum of two overlapping Gaussians, (Fig. 

1), 

fix) = 1.03e-2-2^(^+°-^6)' + 1.03e-2-25(-o.86)^ (7) 

Note that we determined 60 values of f{x) for each Cn, so that the curve in 
Fig. 1 contains 180 points. When x ^ 1, /(x) can be appoximated by 

fix) ~ e-'-''^' (8) 

Therefore, since is proportional to f{x), ns exhibits a stretched exponential 
decay for large S: ns ~ e~^^^'' , in contrast to the exponential decay reported 
in the literature'-^'^-' based on the solution of the percolation problem on the 
Bethe lattice. 

In order to calculate h{x), we first computed the product of h{x)f{x). 
From Eq. (|l|) and Eq. @, we have 

h{Cn - CnJfiCn - C^J = ns^^^ (9) 



The right hand side of Eq. (H) is just the average of Rs/ S'^^^^'^f for fixed S, 
which can be easily calculated numerically. Our numerical results show then 
that h{x)f{x) is proportional to f{x), Eq. (0), (Fig. 2), 

fix)hix) = 0.39e-2-25(-+0-86)^ + Q_3gg-2.25(.-0.86)2 (^q) 

SO that h{x) is essentially a constant, 0.37. 

We have also obtained the first moment of the trajectory length distribution 
as a function of Cr—Cr^. Our numerical calculations show the following power 
law behavior, 

<s>={CR-CRy^ (11) 

where 7 = 2.00 ± 0.01, (Fig. 3), in good agreement with the exact result 
7 = (3 — r)/(T = 2. Our numerical results for the mean square displacement 
of the trajectories also show a power law behavior, 

< >~ {Cr - CrX^ (12) 

where p = 2.33 ± 0.01, (Fig. 4), in good agreement with the exact result 
p = (2 + 2/df - t)Io = 1/a = 7/3 = 2.333. 

As a check, we also used a different method, the so called histogram 
method, to calculate the scaling function f{x) for different concentrations 
from data computed at a given concentration Cr. A similar method has been 
used by Leath*^^'^°^ to study the scaling behavior of percolation clusters and by 
Ferrenberg and Swendsen^^'^^ to study the critical behavior of the Ising model. 
This calculation is based on the following observation. Each trajectory that is 
generated at one concentration of right rotators (Cr) and left rotators (Cl), 
can also be generated at other concentrations. However, the probability with 
which the trajectory is generated depends on the concentrations, 

P{Nr,Nl)=AC^''C^^ (13) 
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where A is the normahzation factor for P{Nji, Nl) and and Nl are the 
number of right rotators and left rotators contained in the trajectory, respec- 
tively. A can also be defined as the total number of all different closed trajecto- 
ries generated on the lattice, so that A does not depend on the concentration. 
The probability to generate the same trajectory at another concentration 
and is 

= PiNn,N,)^{^)^H^)^^ (14) 

where A' is the normalization factor for P'{Nji, Nl). Note that if C'^/Cr > 1, 
then C'l/Cl < 1, or vice versa. To evaluate (CyCL)^^ and (C'l/Cl)^'^ 
separately for large Nji and is not feasible on the computer, since one 
term is too large and the other term is too small. In our simulations, their 
product is computed from an alternative expression, e^^ MC'ii/Cr)+Nl in(cycL)^ 
which works quite well. Here, we emphasize that, although in theory A and 
A' should be the same, for numerical simulations, A and A' are not equal due 
to systematic errors in the Monte Carlo samplings, so that they have be to be 
determined from Eq. (|T^) and Eq. ([T^, respectively, ns a.t concentration 
can be obtained from P'{Nji, Nl) as, 

^5 = E E P'i^R, NL)SiS - S'iNR, Nl)) (15) 

where S'^Nr, Nl) is the length of a trajectory containing Nr left rotators and 
Nl left rotators. Once n5 is obtained, f{x) can be calculated for any other 
concentration according Eq. (j^). 

The standard data that we used were computed at the concentration Cr = 
0.51. From these we generated the scaling functions f{x) for C'^ = 0.52, 0.53, 
0.54, 0.55, all of which fall on the same curve (Fig. 5). Thus this method is 
consistent with the first one. These results also show that the critical region 
in concentration is at least as large as 0.05. 
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2.2 The partially occupied lattice 

For C < 1, the particle trajectories cannot be mapped onto a percolation 
problem anymore, because the trajectories can cross themselves. It was found 
by Cohen and Wang^^^^ that there are two nonlinear critical lines symmetric 
with respect to the line Cr = Cl, merging at Cr = Cl = 1/2, (Fig. 6). 
Moreover, we found that they appear to be tangent to the line C = 1. Previous 
numerical results have shown that the trajectory fractal dimension df and the 
exponent r along the critical lines are the same as that in the case Cr = Cl = 
1/2. Here we show numerically, however, that in addition to the exponent a 
which we have investigated in the previous section, a different exponent a' also 
appears when C < 1, suggesting a new scaling behavior. 

For C < 1, the simplest way to study the critical behavior near criticality 
is to approach the critical point Cr^ = Cl^ = 1/2 along the line Cr = Cl, 
which is perpendicular to the critical line at Cr = Cl = 1/2. We assume that 
the cluster size distribution has a form similar to Eq. ([^), 

ns = S-^+'fiil-C)S'^') (16) 

where f'{x) is a new scaling function in contrast to f{x), which we have 
obtained before and x is defined as (1 — C)S°'' . Since x cannot be negative, 
f'{x) must be asymmetric with respect to x = 0, in contrast to f{x) of Eq. (|^. 
The exponent r, however, must be the same as before, since for C = 1, the 
previously established relation ra^ ~ 5^^+^ must be recovered. 

The computation of the new scaling function f'{x) is much more time con- 
suming than that for the rotator model on the fully occupied lattice, although 
the procedure is similar. The reason is the following. When C is small, the 
moving particle becomes less frequently scattered, therefore, on average, it 
takes longer for the particle trajectory to get closed. On the other hand, if 
C is close to one, i.e. near criticality, the particle trajectories also close very 
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slowly. On the line Cr = Cl, we found that the trajectories close fastest at 
Cr = Cl = 0.45, when essentially, all trajectories are closed after 2^^ time 
steps (Fig. 7). However, this concentration is not useful to calculate the 
scaling function, because it is not in the critical region. 

As we explained in the introduction, if the critical behavior appears only 
for very large trajectories, direct calculation of 7' and p' — the exponents 
corresponding to 7 and p, respectively — from equations Eq. (^) and Eq. (^) 
becomes very difficult, because the large trajectories are too large to generate 
on the computer. On the other hand, we can still compute the scaling function 
f'{x) by making a cutoff, i.e. trajectories whose length is bigger than a fixed 
number are disregarded. In our simulations, we took this fixed number or 
cutoff as sixteen million time steps. By choosing a' = 1.6, we find that the 
scaling functions obtained for different concentrations, C = 0.99, 0.985, 0.98, 
all collapse to a single curve, (Fig. 8), 

fix) = 0.475e-(^-^^^^°"'") (17) 

where for each C, 25 points have been computed, so that the curve in Fig. 8 
contains 75 points. Note that since x ~ S^-^, f\x) decays as a super-exponential 
function of S", ng ~ e"'^^ '^. In fact, eq.(17) obtains for values of a' = 1.6 ±0.3. 

The other scaling function h'{x) is obtained from the calculation of f'{x)h'{x), 
which we found again to be proportional to f'{x), (Fig. 9), 

/'(x)/i'(x) =0.115e-(^-^^^^°~'") (18) 

so that h'{x) = 0.115/0.475 = 0.24 independent of x. The exponents 7' and p' 
can be obtained from 7' = (3 - r) /a' = 0.54 and p' = (2 + 2/dj - r) /a' = 0.63, 
respectively, both much smaller than those for C = 1. 
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3 Critical region for the mirror model on the 
square lattice 



It has been shown before, by mapping the mirror model to a percolation prob- 
lem, that C = 1 is a critical lme^^'^\ i.e. the size distribution and the fractal 
dimension of the trajectories have the same power law behavior as for the cor- 
responding percolation problem at criticality, ns ~ S~'^'^^ ~ 5""^/^ or r = 15/7 
and df — 7/4, respectively. However, the trajectory size average < S > and 
the mean square displacement < > are both divergent along the critical 
line. Therefore, the exponent a does not appear here. Although the mirror 
model can be mapped onto a quasi-rotator model which exhibits the same 
scahng behavior as that for the ordinary rotator model at criticality^^^ , the 
quasi-rotator model cannot be used to calculate a, 7 and p, because at present 
we do not know how to generate the particle trajectories directly rather than 
deducing them from the mirror model. 

For C < 1, the mirror model cannot be mapped onto a percolation problem 
anymore. It was found numerically that the critical exponents are drastically 
changed from those at C — 1; ns ~ (InS*)"^ and df — 2 with logarithmic 
corrections'^^'®-'. This new critical behavior exists in the whole {Cr, Cl) plane, 
where super-diffusion occurs^^\ except at the three boundary lines C = 1, 
Cr — and C^, = 0, where the moving particle simply zig-zags to infinity. 
Since < 5" > and <R^> are both divergent everywhere, a cannot be defined 
either. 

4 Critical region for the rotator and the mir- 
ror model on the triangular lattice 
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4.1 The fully occupied lattice 

Since the rotator model and mirror model are equivalent on the triangular 
lattice*^^'^'^) , we only study the rotator model here. Similar to the rotator 
model on the square lattice, the rotator model on the triangular lattice can 

be mapped onto a percolation problem*^^'^'''^ a site percolation problem here, 
for which r — 15/7, df — 7/4 and cr — 3/7, while the critical concentration is 
Cr^ = Cl^ = 1/2. 

Our numerical calculations of the scaling function f{x) were carried out at 
Cr = 0.47, 0.48 and 0.49, respectively. The scahng functions for these three 
concentrations collapsed to a single curve, which could be fitted again to a 
double Gaussian, (Fig. 10), 

f{x) = o.94e-i-«°(^+°-^^)' + o.94e-^-«°(^-°-^6)' (19) 

where for each Cr, more than 60 points were computed and the curve comprises 
therefore more than 180 points. As before, f{x) must be symmetric with 
respect to x = 0, the reason being that the probability to generate the same 
trajectory is invariant under the transformation of interchanging Cr and C^. 
Eq.(7) and Eq.(19) are very similar escept that the constants are different. 
When a; ^ 1, /(x) ~ e~^^, so ~ e"'^^^^ which is the same for the square 
lattice, i.e. the number of trajectories decays with a stretched exponential law 
in S. 

The calculation of the scahng function h{x) is similar to that for the square 
lattice. We found numerically, f{x)h{x) is proportional to f{x), (Fig. 11), 

f{x)h{x) = o.28e-^-«°(^+°-^^)' + o.28e-^-«°(^-'^-^^)' (20) 

so that h{x) is a constant, 0.28/0.94 = 0.30. 

The average trajectory size diverges as one approaches criticality with an 
exponent j, < S >~ {Cr — CrJ'"^. Our numerical simulations show that 
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7 = 2.0 ± 0.01, (Fig. 12), in good agreement with the exact result, 7 = 2. 
The mean square displacement of the particle trajectories diverges with an 
exponent p, < Ft' >~ {Cr — Cr^)~'' , the value of p obtained from our numerical 
simulations is 2.33 ± 0.01 (Fig. 13), also in good agreement with the exact 
result 7 = 7/3. 

Using the histogram method, we computed the scaling functions also from 
the standard data obtained at the concentration Cr = 0.51, for the other 
concentrations Cr = 0.52, 0.53, 0.54, 0.55. All the scaling functions fall on 
the same curve, Eq. (|19|), (Fig. 14), showing that this method is consistent 
with the first one. 

4.2 The partially occupied lattice 

For C < 1, the rotator model and the mirror model can be still mapped 
onto each other. However neither model can be mapped onto a percolation 
problem now. Here we only consider the rotator model. It was found before 
that there is only one critical line Cr = Cl rather than two critical lines as in 
the case of the rotator model on the partially occupied square lattice^^'^). The 
critical exponents r and df along the critical line are found numerically to be 
the same as those at the critical point, Cr = Cl = 1/2, viz. 15/7 and 7/4, 
respectively*^^'^' . 

To study the critical behavior near criticality, we first need to choose the 
direction in which we approach criticality. The most obvious choice is the 
direction perpendicular to the critical line. Although the rotator model cannot 
be mapped onto a percolation problem, we found that, if we choose cr = 3/7, 
the scaling functions computed at different Cr collapse very well onto a single 
curve. This suggests that randomly distributed empty sites on the lattice are 
irrelevant for the critical behavior. The scaling functions were computed in the 
critical region near Cr^ = Cl, = 0.425 at three different Cr, Cr = Cr, + 0.01, 
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Cn = Cr^ + 0.015, Cr = Cr^ + 0.02 along the line C = 0.85. We found that 
these scaling functions could be fitted to a double Gaussian, (Fig. 15), 

fix) = 1.05e-i-^«(^+^-i9)' + 1.05e-i-i«(^-i-i9)' (21) 

where we computed 80 points for each Cr, leading to 240 points on the curve. 
From Eq. (|l^)and Eq. (pT]), one can see that all three constants contained in 
/(x) are concentration dependent. 

The product of f{x)h{x) is again found to be proportional to /(x), (Fig. 
16), 

f{x)h{x) = 0.30e-i-^«(^+^-i9)' + 0.30e-i-^«(^-^-i9)' (22) 

so that h{x) = 0.30/1.05 = 0.28. 

To investigate whether the critical exponent a has the same value along 
the critical line, we also computed < S > and < i?^ > in the critical region 
around Cr^ = Cl^ = 0.45 and Cr^ = Cl^ = 0.4, respectively. In both cases, 
we found, from our numerical simulations, that the exponents 7 = 2.00 ± 0.02 
and p = 2.33 ± 0.02 are very close to their values at Cr^ = Cl^ = 0.5: 2 and 
7/3, respectively, (Fig. 17-18), suggesting that the critical behavior along the 
critical line Cr = Cl belongs indeed to the same universahty class. 

5 Conclusion 

We end with a few remarks. 

1. We have studied the critical behavior of particle trajectories on both the 
square lattice and the triangular lattice in the critical region, for both C = 1 
and C < 1. Our methods were based on calculations of two scaling functions 
f{x) (which yields a) and h{x), as well as the two exponents 7 and p from 
< S> and <R^> respectively. From these two scaling functions, one derived 
and verified numerically the two scaling relations, 7 = (3 — r)/cr and p = 
(2 + 2/df — t)/(t. Our results are summarized in table I. 
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2. For C = 1, the rotator model on the square lattice, and the rotator and the 
mirror model on the triangular lattice can all be mapped onto a percolation 
problem. Our numerical results show that near Cr^ — Clq — 1/2, f{x) can 
be fitted to a sum of two overlapping Gaussians, i.e. a double Gaussian. Our 
numerical calculations of < 5" > and < > show that 7 = 2.00 ± 0.01 and 
p = 2.33 ± 0.01, respectively , in excellent agreement with the exact results. 
Weinrib and Trugman*^^^) have argued that the scaling function f(x) for the 
perimeter of percolation clusters should be proportional to the scaling function 
g{x) for the percolation clusters themselves. It would be interesting to check 
this for the f{x) we obained here. 

3. For C < 1, however, the critical behavior of the rotator model on the 
square lattice and that of the rotator and mirror model on the triangular lattice 
are very different. None of these models can be mapped onto a percolation 
problem. For the square lattice, we obtained numerically a new exponent 
a' — 1.6 ± 0.3 which is significantly larger than the usual a — 3/7. This 
indicates that the typical size of closed trajectories upon approach to criticality 
increases much slower to infinity along the line Cr = Cl than along the line 
Cr + Cl — 1. The new scaling function f'{x), corresponding to the exponent 
a' — 1.6, appears to be an exponential function, rather than a double Gaussian 
as for C = 1, but the origin of this exceptional critical behavior is unclear to 
us. For the triangular lattice, however, we found numerically, that a is very 
close to 3/7 which also obtains for C = 1, and that the scaling function 
f{x) can still be described by a double Gaussian. This suggests that the 
critical behavior in the direction perpendicular to the critical line Cr — Cl 
is in the same universality class. For dilute thermal systems, it is also found 
that the critical behavior is independent of the density of impurities if the 
specific heat exponent a is negative*^^^\ but to what extent this is similar to 
the independence of the number of empty sites found here, is unclear to us. 
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4. The scaling function h{x) appears to behave hke a constant in all cases, 
indicating that there are no noticeable corrections to the universal fractal di- 
mension df = 7/4. This simply means that the gyration radius of the extended 
trajectories is independent of the degree of occupation of the lattice by scat- 
terers and only trivially dependent on the nature of the lattice. 

5. We have studied in the previous paper the scaling behavior of extended 
particle trajectories as their sizes approach infinity at criticality. For exam- 
ple, we have obtained for the rotator model on the fully occupied square 
lattice that < S/N — 3/2 >~ N"^-^"^, where the exponent is very close to 
1 — a — 4/7 = 0.571, with a characterizing the trajectory size distribution 
in the critical region near criticality. We surmise that this may not be a co- 
incidence and conjecture that there is a possible relation between the scaling 
behavior of extended closed trajectories at criticality and in the critical region 
near criticality. If so, the question arises whether also other exponents such 
as 7 and p could be determined from the scaling behavior at criticality, rather 
than as we have done here, in the critical region. 

6. Finally, we remark that the diffusive behavior of the moving particle away 
from criticality, which has been studied extensively before^^''^^ is in fact con- 
trolled by the scaling function f{x). The time dependent diffusion constant 
D{t) is defined as D{t) =<i?^(i) > /{At) = E R%ns/At, where the summation 
is taken over all closed trajectories upto length t. For [Cr — C^Jf^ <^ 1, when 
f{x) is essentially a constant, D{t) is then virtually the same as at criticality, 
while for {CR—CR^f^ ^ 1, when f{x) is very small, D{t) decays to zero asl/t. 
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Figure captions 



Fig. 1. The scaling function f{x) vs x for the rotator model on a fully 
occupied square lattice, computed at C = Cr^ + 0.01 (o), C = Cr^ + 0.02 
(+) and C = Cr^ + 0.03 (□). The curve is described by a double Gaussian, 
f{x) = 1.03e-2-25(a:+o.86)2^i Q3g-2.25(x-o.86)2 The deviations between f{x) and 

the numerical data near a; = in this and later similar figures are due to the 
failure of scaling for small trajectories. 

Fig. 2. The scahng function f{x)h{x) vs x for the rotator model on a fully 
occupied square lattice, computed at C = Cr^ + 0.01 (o), C = Cr^ + 0.02 

(+) and C = Cr^ + 0.03 (□). The curve is described by a double Gaussian , 

f{x)h{x) = 0.39e-2-25{-+o.86)^ +o.39e-2-25(-o.86)2^ 

Fig. 3. In < S'> vs — ln(Cii — Crc) ^r the rotator model on the fully occupied 
square lattice. The slope of the fitting hne is 2.00. 

Fig. 4. In < i?^ > vs — ln(Cij — Cr^) for the rotator model on the fully occupied 
square lattice. The slope of the fitting hne is 2.33. 

Fig. 5. The scaling function f{x) vs x for the rotator model on the fully 
occupied square lattice, obtained by using the histogram method. Computed 
for C ^ Cr^ + 0.02 (o), C = Cr^ + 0.03 (+), C = Cr^ + 0.04 (□) and 
C = Cr^ + 0.05 (x). The curve is described by a double Gaussian, f{x) — 

Q3g-2.25(a;+0.86)2 _^ Q3g-2.25(a;-0.86)2 

Fig. 6. Part of the phase diagram for the rotator model on the square lattice 
obtained from numerical simulations for 0.65 < C < 1. The lines are drawn 
to guide the eye. 

Fig. 7. ns vs 5", at concentrations Cr^Cl^ 0.5 (❖), Cr^Cl^ 0.485 (+), 
Cr = Cl = 0.45 (□) and Cr = Cl = 0.40 (x). 

Fig 8. The scaling function f'{x) vs x for the rotator model on the partially 
occupied square lattice, computed at C = Ci?^ + 0.02 {<>)^C — Cij^ + 0.025 (+) 
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and C = Cr^ + 0.01 (□). The curve is described by an exponential function, 
f'{x) = 0.475e-°-^65^i'^"'^ 

Fig 9. The scahng function f'[x)h'{x) vs x for the rotator and mirror model 
on the partially occupied square lattice, computed at C = Cr^ + 0.02 (o), 
C = Cij^ + 0.025 (+) and C = Cr, + 0.01 (□). The curve is described by an 
exponential function, f'{x)h'{x) = O.llSe"^--'^^^^-'^'^ *^. 

Fig 10. The scaling function f[x) vs x for the rotator and mirror model 
on the fully occupied triangular lattice, computed at C = Cr^ + 0.01 (o), 
C ^ Cr^ + 0.02 (+) and C ^ Cr^ + 0.03 (□). The curve is described by a 
double Gaussian, f{x) = 0.94e-i-«(^+o-96)' + 0.94e-i-«(^-°-96)'. 
Fig 11. The scaling function f{x)h{x) vs x for the rotator and mirror model 
on the fully occupied triangular lattice, computed at C = Cr^ + 0.01 (o), 
C ^ Cr^ + 0.02 (+) and C ^ Cr^ + 0.03 (□). The curve is described by a 
double Gaussian, f{x)h{x) = 0.28e-i-«(^+°-96)' + 0.28e-i-«(^-°-96)'. 
Fig. 12. In < 5" > vs — In |Cij — CrJ\ for the rotator and mirror model on the 
fully occupied triangular lattice. The slope of the fitting hne is 2.00. 
Fig. 13. In < i?^ > vs — In \Cr — CrJ\ for the rotator and mirror model on 
the fully occupied triangular lattice. The slope of the fitting line is 2.33. 
Fig. 14. The scahng function f{x) vs x, for the rotator and mirorr model on 
the fully occupied triangular lattice , obtained by using the histogram method, 
computed for C = Cr^ + 0.02 (o), C = Cr^ + 0.03 (+), C = Cr^ + 0.04 (□) 
and C = Cr^ + 0.05 (x). The curve is described by a double Gaussian, 
f{x) = o.94e"^-^^(^+°-^^)'' + o.94e"^•^^(^"°•^'^)^ 

Fig 15. The scaling function f{x) v.s x for the rotator and mirror model on 
the partially occupied triangular lattice, computed at C = Cr^ + 0.01 (o), 
C = Cr^ + 0.015 (+) and C = Cr^ + 0.02 (□). The curve is described by a 
double Gaussian f{x) = i.05e-i-i^(^+^-i9)' + l.Obe-^-^^^''-^-^^^' 
Fig 16. The scaling function f{x)h{x) vs x for the rotator and mirror model 
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on the partially occupied triangular lattice, computed at C = Cr^ + 0.01 (o), 
C = Cr^+ 0.015 (+) and C = Cr^ + 0.02 (□). The curve is described by a 
double Gaussian f{x)h{x) - O.SOe-^-^^^^+i-^^)' + O.SOe-^-^^^^-^-^^)' 
Fig. 17. In < 5" > vs — \n.\CR — CrJ\ for the partially occupied triangular 
lattice, for C = 0.9 (❖)and C = 0.8 (+). The slope of the fitting hues is 2.00. 
Fig. 18. In < i?^ > vs — ln|CR — CrJ\ for the partially occupied triangular 
lattice, for C = 0.9 (o)and C = 0.8 (+). The slope of the fitting lines is 2.33. 
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Table I: Scaling functions and critical exponents 





Lattice 






Square 


Triangular 


Scaling Functions 




C = 1 


C < 1 


C = 1 


C < 1 


f{x) (double Gaussian) 


A, 

«i 
fli 


1.03 
2.25 
0.86 




0.94 
1.80 
0.96 


1.05 
1.18 

1.19 

(C = 0.85) 


h{x) (constant) 




0.38 




0.30 


0.28 


f'{x) (exponential) 


A', 

'A 




0.475 
1.65 -10-^ 








h'{x) (constant) 






0.24 






Exponents 












r 

df 
a {a') 
7(7') 

P ip') 




15/7 
7/4 
3/7 
2 

7/3 


15/7 

7/4 

1.6 

0.54 

0.63 


15/7 
7/4 
3/7 
2 

7/3 


15/7 
7/4 
3/7 
2 

7/3 



Scaling functions and critical exponents obtained for closed trajectories on the 
square and triangular lattice. f{x) = Ai exp[— ai(a; — ai)^] + Ai exp[— ai(a; + 
aiY] and f'{x) — A[ exp[ —a[{x — a'l)]. The critical exponents for C = 1 are 
known exactly, those for C < 1 only numerically. The square lattice for C < 1 
behaves exceptional if the critical point Cr^ = Cl^ = 1/2 is approached along 
the line Cr = Cl, rather than along the line C = 1; the primed quantities 
refer to this case. 
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